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The Atlas Mountains, not the Strait of Gibraltar, 
as a biogeographic barrier for Mauremys leprosa 
(Reptilia: Testudines) 

UwE Fritz, Guido Fritzsch, Edgar Lehr, Jean-Marc Ducotterd & Anke Muller 

Abstract. Sequence data of the mitochondrial cytochrome b gene of stripe-necked terrapins {Mauremys 
leprosa) have heen compared from localities north and south of the Atlas Mts. (Ceuta; Morocco) and from 
Donana National Park (Spain). A low maximum sequence divergence (approximately 1%) corresponds 
to two major clades; one is represented hy localities to the north of the Atlas Mts. and in Donana National 
Park and the other by localities to the south of the Atlas Mts. Differentiation between populations north 
and south of the Atlas Mts. is much more pronounced than that found between samples from each side 
of the Strait of Gibraltar. These findings suggest that the Strait of Gibraltar is, in contrast to the Atlas Mts., 
not a significant barrier to gene flow in stripe-necked terrapins. The major clades could reflect taxonomic 
segregation between populations north and south of the Atlas Mts. Sequences from Marrakech (corre¬ 
sponding to M. 1. rnarokkensis), Ceuta, and the Donana National Park (M. 1. leprosa) are only weakly 
differentiated. South of the Atlas Mts. we found no consistent differences between samples from catchment 
basins of the Oued Draa (M. 1. vanmeerhaghei) and the Oued Noun {M. 1. saharica). Our findings imply 
that taxonomic differentiation within M. leprosa is currently overestimated. 
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Introduction 

The Strait of Gibraltar has been highlighted 
as an important biogeographic barrier that 
separates North African and Iberian herpeto- 
faunas (Busack 1986). Recent mtDNA se¬ 
quence data corroborate the barrier status of 
the Strait of Gibraltar for some amphibians 
and reptiles (Alytes; Fromhage et al. 2004; 
Discoglossus; GarcIa-ParIs & Jockusch 1999, 
Fromhage et al. 2004; Pelobates; Garc1a-Pa- 
rIs et al. 2003; Acanthodactylus: Harris et al. 
2004), but not for others (Pleurodeles waltl; 
Batista et al. 2004, Carranza & Arnold 2004, 
Veith et al. 2004; Podarcis: Harris et al. 
2002). The stripe-necked terrapin {Maure¬ 
mys leprosa), widely distributed on both 
sides of the Strait of Gibraltar, was long 
thought to be monotypic (Busack & Ernst 
1980). ScHLEicH (1996) and Bour & Maran 
(1999), however, described seven subspecies 
from Morocco, suggesting considerable vari¬ 


ation. Recent molecular investigations pro¬ 
posed M. leprosa as sister taxon to other 
West Palearctic Mauremys species (M. caspi- 
ca, M. rivulata) and to several East Asiatic 
taxa (Barth et al. 2004, Spinks et al. 2004). 
The split between M. leprosa and the other 
taxa possibly dates back to the Late Oligo- 
cene or Early Miocene (Barth et al. 2004), 
approximately matching the age of the old¬ 
est known West Palearctic Mauremys fossils 
(Lapparent de Broin & van Diik 1999, Lappa- 
RENT DE Broin 2001, Her vet 2004). Recent 
subspecies descriptions, the long independ¬ 
ent evolutionary history of M. leprosa, and 
the presence of old Mauremys-\\ke fossils in 
West Europe and North Africa (Lapparent de 
Broin & van Diik 1999, Lapparent de Broin 
2000, 2001, Hervet 2004) suggest that popu¬ 
lations on each side of the Strait of Gibraltar 
could be old and well-differentiated. We tes¬ 
ted this hypothesis by comparing sequences 
of the mitochondrial cytochrome b gene 
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Fig. 1. Geographic origin of samples: (1) Donana 
National Park; (2) Ceuta; (3) SE Marrakech; (4) 
Oued Dades, NE Ouarzazate; (5) Oued Draa, SE 
Ouarzazate; (6) Oued Tata, S Tata; (7) Oued Noun, 
E Guelmine; (8) Oued Noun, N Tiliouine. For exact 
localities see Table 1. 


(cyt b) of M. leprosa from southwest Spain 
(Donana National Park), the Spanish exclave 
Ceuta (North Africa), and Morocco (locali¬ 
ties north and south of the Atlas Mts). Cyt b 
is known to he phylogenetically informative 
in testudinoid chelonians (Caccone et al. 
1999, Lenk et al. 1999, Feldman & Parham 

2002, Palkovacs et al. 2002, Stephens & Wiens 

2003, Barth et al. 2004, Spinks et al. 2004). 


Materials and Methods 

Sampling 

We studied \2> Mauremys leprosa from Spain 
and Morocco (Fig. 1, Tah. 1). One specimen 
was included in a previous study (Barth et al. 
2004) and most were captured by one of us 
(J.-M. D.). Individuals were released after 
taking blood samples, measurements, and 
photographs. Morphological characters and 
habitats are described in Ducotterd & Bour 


(2002). Blood samples were preserved in 
98% pure ethanol or EDTA buffer (0.1 M 
Tris, pH 7.4, 10% EDTA, 1% NaF, 0.1% 
thymol) and stored at -80 °C. Voucher pho¬ 
tographs, blood and remaining DNA samples 
are permanently housed in the tissue collec¬ 
tion of the Museum of Zoology Dresden 
(Museum fiir Tierkunde, MTD T). 


Sequencing 

Total genomic DNA was extracted from 
small quantities of blood following Gustin- 
ciCH et al. (1991). For amplification of cyt b 
we used the primers mt-a-neu (5’-CTC CCA 
GCC CCA TCC AAC ATC TCA GCA TGA 
TGA AAC TTC G-3’) of Lenk & Wink (1997) 
and H-15909 (5’-AGG GTG GAG TCT TCA 
GTT TTT GGT TTA CAA GAC CAA TG-3’) 
of Lenk et al. (1999). PCR was conducted in 
an Eppendorf Mastercycler. Cycling proce¬ 
dure was: initial denaturation 5 min at 95 °C; 
40 cycles, denaturation 1 min at 95 °C, pri¬ 
mer annealing for 1 min at 50 °C, extension 
for 2 min at 72 °C, and stop reaction for 10 
min at 72 °C. PCR products were purified 
using the Invisorb Spin PCRapid Kit from 
Invitek. Sequencing primers were Mau-F (5’- 
CTA GGC CTC ATC TTA ATA CT-3’) and 
Ri-neu (5’-GTG AAG TTG TCT GGG TCT 
CCT AG-3’). Sequencing was performed on 
an ABI3100 Genetic Analyzer (Hitachi) with 
the following conditions for 25 cycles: 10 s 
at 96 °C; 5 s at 50 °C; and 4 min at 60 °C. 
Obtained sequences were approximately 
1,100 bp long. 


Phylogenetic analysis 

MEGA 2.1 (Kumar et al. 2001) was used to 
estimate genetic distances and calculate se¬ 
quence statistics. Gene fragments were 
aligned with ClustalX v. 1.83 (Thompson et al. 
1997) and default parameters. Aligned se¬ 
quences (including gaps) comprised 993 bp. 
Mauremys rivulata, retrieved from GenBank 
(AJ564455; Barth et al. 2004), served as the 
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Fig. 2. Neighbour-Joining tree of the cyt b dataset for Mauremys leprosa, rooted with M. rivulata. Clusters 
1 and 2 represent localities north of the Atlas Mts. and from both sides of the Strait of Gibraltar, clusters 
3 and 4 localities south of the Atlas Mts. Large numbers at the leaves are MTD T numbers (Tab. 1); numbers 
in black circles denote collection sites (Fig. 1). The sample from the Oued Draa catchment area with the 
haplotype occurring otherwise in the Oued Noun region is asterisked. Small numbers at the nodes are, 
from left to right: bootstrap values of 10,000 trees for NJ, bootstrap values for ML (TBR swapping 
algorithm), posterior probabilities for MrBayes, and bootstrap values for MP (1,000 resamplings). Branch 
lengths correspond to number of nucleotide changes. 


outgroup. Different procedures were used for 
estimating phylogeny. To find the most ap¬ 
propriate model for DNA nucleotide substi¬ 
tution, we performed a hierarchical likeli¬ 
hood ratio test with Modeltest v. 3.5 (Posada 
& Crandall 1998). For the dataset, HKY 
(Hasegawa et al. 1985), a TRatio = 7.2040, 
and equal rates provided the best fitting of 56 
tested models. These parameters were used to 
perform Maximum Likelihood (ML) and 
Bayesian analysis. The ML analysis was cal¬ 
culated with PAUP* 4.0b 10 (Swofford 2002) 
using the heuristic search method with 10 
random stepwise additions and TBR branch 
swapping for 1,000 bootstrap replicates. 
Bayesian analyses were performed with 
MrBayes V3.0b4 (Huelsenbeck & Ronquist 
2001), which was used to run 1,000,000 


generations, with a sampling frequency of 10 
generations. From the 100,000 trees found, 
the first 1,000 were discarded. Neighbour- 
Joining (NJ) trees (Saitou & Nei 1987) were 
constructed with MEGA 2.1 and PAUP*. We 
chose the model of Kimura-2-parameter (Ki- 
MURA 1980), a pairwise deletion and 10,000 
bootstrap replicates. Maximum Parsimony 
(MP) analyses were performed with PAUP* 
using the heuristic search method with 10 
random stepwise additions and the TBR 
branch swapping option. Bootstrap analyses 
(Felsenstein 1985) were used to examine ro¬ 
bustness of resulting bifurcations within the 
trees. The MP tree was tested with 1,000 
replicates. 

In addition, we used Bio-Neighbour-Joi- 
ning (Gascuel 1997) and Neighbour-Net 
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Fig. 3. Neighbour-Net of the cyt b dataset foi Mauremys leprosa, rooted with M. rivulata. Small numbers 
along branches denote bootstrap values (1,000 replications), indicating probability for respective pathway. 
Large numbers represent MTD T numbers (Tab. 1). Branch length is proportional to number of divergent 
characters supporting this pathway. For further explanation, see Figure 2. 


Locality 

n 

Referred 

subspecies 

MTD T 

EMBL accession 
number 

1 Spain: Donana National Park, 

Laguna Dulce 

2 

leprosa 

2506,3305 

AJ564457, AJ877040 

2 Spain: Ceuta 

1 

7 

1299 

AJ877039 

3 Morocco: SE Marrakech: near Ait-Ourir, 
31°27.891N, 7°46.724W; 672 m 

2 

marokkensis 

164-165 

AJ877025-AJ877026 

4 Morocco: NE Ouarzazate: Sidi-Flah, Oued 
Dades, 31°00.688N, 6°29.931W; 1158 m 
(MTD T 763, 768) and 31°00.695N, 
6°30.808W; 1177 m (MTD T 769) 

3 

vanmeerhaghei 

763,768-769 

AJ877024, 

AJ877027-AJ877028 

5 Morocco: SE Ouarzazate: Tamnougalt, 

Oued Draa, 30M0.183, 6°22.881W, 917 m 

2 

vanmeerhaghei 

770-771 

AJ877029-AJ877030 

6 Morocco: 18 km S Tata: El-Khemis, 

Oued Tata, 29°35.506, 8°00.009W, 494 m 

2 

vanmeerhaghei 

774, 779 

AJ877031-AJ877032 

7 Morocco: 17 km E Guelmine: Oued 

Noun, 28°58.433N, 9°54.231W, 300 m 

2 

saharica 

780, 784 

AJ877033-AJ877034 

8 Morocco: N Tiliouine: Oued Noun 

canyon, 29°05.115N, 10°15.140W, 122 m 

4 

saharica 

785-788 

AJ877035-AJ877038 


Tab. 1. Mauremys leprosa samples, localities and EMBL accession numbers. MTD T: tissue collection. 
Museum of Zoology Dresden (Museum fiir Tierkunde). 
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Cluster 1 

Cluster 2 

Cluster 3 

Cluster 4 

Mauremys rivulata 

Cluster 1 

- 

3 

10 

9 

75 

Cluster 2 

0.3 

- 

9 

8 

72 

Cluster 3 

1.0 

0.9 

- 

5 

73 

Cluster 4 

0.9 

0.8 

0.5 

- 

76 

Mauremys rivulata 

7.6 

73 

7.4 

7.7 

- 


Tab. 2. Divergence (p-distance) between clusters of cyt b sequences (993 bp) of Mauremys leprosa and 
the outgroup M. rivulata. Number of changes above diagonal, sequence divergence (%) below diagonal. 
Average values are given for cluster 1. Cluster 1: Ceuta and Donana National Park; cluster 2: Marrakech; 
cluster 3: catchment area of Oued Draa; cluster 4: catchment areas of Oued Noun and Oued Draa. 


(Bryant & Moulton 2004) as implemented in 
the program SplitsTree v. 4betal0 (Huson 
1998, Huson & Bryant 2004). Neighbour-Net 
is a powerful tool for testing the consistency 
of phylogenetic hypotheses (Bryant & 
Moulton 2004). 


Results 

Nine hundred ninety-three base pairs of the 
Mauremys leprosa cytochrome b gene were 
unambiguously aligned for phylogenetic 
purposes. Within this alignment, 82 posi¬ 
tions were variable and 13 were parsimony 
informative. These differences correspond to 
five weakly-differentiated haplotypes. The 
greatest sequence divergence of approxi¬ 
mately 1% (p-distance) exists between lo¬ 
calities north of the Atlas Mts. plus south¬ 
west Spain versus localities south of the 
Atlas Mts. (Tab. 2). 

All phylogenetic methods resulted in the 
same topology. The Neighbour-Joining (NJ) 
tree (Fig. 2) and the Neighbour-Net (Fig. 3) 
serve as examples. In the NJ tree bootstrap 
and probability values for Maximum Likeli¬ 
hood, Bayesian and Maximum Parsimony 
analysis are included. In all presentations, M. 
leprosa sequences are monophyletically 
rooted with M. rivulata. A major split, sup¬ 
ported by moderate to high bootstrap and 
probability values, occurs between samples 
from localities north of the Atlas Mts. plus 
southwest Spain versus localities south of 
the Atlas Mts. 


Within the group containing the samples 
from the area north of the Atlas Mts. plus 
southwest Spain, two clades are found. These 
groupings correspond to Ceuta plus Donana 
National Park (MTD T 1299, 2506, 3305) 
and two samples from the area around 
Marrakech (MTD T 764-765), respectively. 
In our sample from Ceuta we found a haplo- 
type identical to that in one individual from 
Donana National Park. Turtles from the area 
south of the Atlas Mts. make up the other 
major branch, consisting also of two clades. 
These clades correspond to individuals from 
the catchment areas of the Oued Noun or the 
Oued Draa, with one exception: one sample 
(MTD T 769) from the Oued Dades, a tribu¬ 
tary of the Oued Draa, shares its haplotype 
with specimens from the Oued Noun. Se¬ 
quence divergence between the Ceuta plus 
Donana clade and the Marrakech clade is 
0.3%, while between the two clades from 
south of the Atlas Mts. the p-distance is 0.5% 
(Tab. 2). 

Because we encountered few phylogene- 
tically informative sites in the cyt b frag¬ 
ment, the branching pattern of the tree-build¬ 
ing methods was tested by Neighbour-Net 
(Bryant & Moulton 2004), a distance-based 
method for constructing phylogenetic net¬ 
works. Unlike a NJ tree, a Neighbour-Net 
(NNet) can represent conflicting phylogene¬ 
tic signal in the data and demonstrate simul¬ 
taneously alternative topologies. NNet al¬ 
lows for reticulations and multifurcations, 
reflecting uncertain domains and alternative 
pathways. If alternative pathways exist or 
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Fig. 4. Mauremys leprosa saharica with pale blue 
eyes (Oued Noun canyon, N Tiliouine, Morocco, 
29°05.115N, 10°15.140W). 


data are ambiguous, branches are represented 
in the graphical network by boxes. In our 
data ambiguities regarding the position of 
the M. leprosa clusters relative to the out¬ 
group M. rivulata exist. The branching pat¬ 
tern of M. leprosa obtained with tree-build¬ 
ing methods is confirmed with NNet however 
(Fig. 3). 


Discussion 

It is generally accepted that North Africa and 
the Iberian peninsula were last in contact 
during the Messinian salinity crisis (Late 
Miocene, approx. 5.5 Mya) when the Medi¬ 
terranean basin dried up due to evaporation 
and the closure of previous Miocene sea 
straits connecting it with the Atlantic. With 
restoration of the connection with the Atlan¬ 
tic, the Mediterranean basin again filled du¬ 
ring a catastrophic event, lasting approxi¬ 
mately 100 years (Hsu et al. 1977, Krijgsman 
et al. 1999, Duggen et al. 2003). Since then. 
North Africa and the Iberian peninsula have 
been separated by shifting sea straits in the 
area of the current Strait of Gibraltar (de Jong 
1998). If North African and Iberian Maur¬ 
emys leprosa populations were separated by 
this vicariant event, genetic differentiation 
should he evident. 

The oldest known northwest African fos¬ 
sils referred to M. leprosa (dated to the Rus- 


cinian, 5.4-3.4 Mya, the Pliocene stage im¬ 
mediately following the Messinian) are from 
Am Boucherit, Constantine province, Alge¬ 
ria. Undetermined Mauremys findings of the 
same age were also excavated in Tunisia 
(Lapparent de Broin 2000). The contemporary 
European fossil species M. gaudryi (Deperet, 
1885), thought to be closely related to M. 
leprosa (Hervet 2004), is known from Perpig¬ 
nan (Department Pyrenees-Orientales, 
France), suggesting that ancestors of extant 
M. leprosa were present in North Africa and 
Europe when the sea strait between the con¬ 
tinents re-opened during the Messinian. 
However, Lapparent de Broin (2001) believes 
that M. leprosa originated in North Africa 
and arrived in Europe not before late Plio¬ 
cene or early Pleistocene. 

We found identical or very similar 
mtDNA sequences on both sides of the Strait 
of Gibraltar, ranging only to 0.3% sequence 
divergence (Tab. 2). Similarly, Mantziou et 
al. (2004) reported little genetic differentia¬ 
tion between populations of M. rivulata oc¬ 
curring on the Greek mainland and Aegean 
islands, and Lenk et al. (1999) reported a 
Moroccan Emys orbicularis cyt b sequence 
resembling haplotypes occurring on the Ibe¬ 
rian peninsula. In each case, rafting or island 
hopping might have prevented development 
of distinct mtDNA lineages. On the other 
hand, we detected pronounced differentia¬ 
tion in M. leprosa occurring south of the 
Atlas Mts. In E. orbicularis the Pyrenees, 
Alps, Dinarides, Carpathians, and the Cauca¬ 
sus Mts. are known to be significant 
biogeographic barriers, separating subspe¬ 
cies that harbour distinct mtDNA lineages 
(Fritz 1992, 1996, 2003, Lenk et al. 1999). 
This suggests that high mountain chains are 
distinctly more important barriers than sea 
straits for freshwater turtles. 

Seven subspecies of M. leprosa that have 
recently been described from Morocco occur 
either north or south of the Atlas Mts. 
ScHLEiCH (1996) described the subspecies M. 
1. atlantica, M. 1. erhardi, and M. 1. werner- 
kaestlei from northern Morocco and M. 1. 
marokkensis from the region of Marrakech. 
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Fig. 5. Juveniles of (a) Mauremys leprosa vanmeerhaghei (Sidi-Flah, Oued Dades NE Ouarzazate, 
31 °00.688N, 6°29.931W) and (b) M. 1. saharica (Oued Noun canyon, N Tiliouine, Morocco, 29°05.115N, 
10°15.140W). Note carapacial patterns. 


From the southern slope and the area south of 
the Anti Atlas, Schleich (1996) described M. 
/. saharica from Oued Noun and M. 1. zizi 
from Oued Ziz, while Bour & Maran (1999) 
described M. 1. vanmeerhaghei from Oued 
Draa. With the exception of M. 1. vanmeer¬ 
haghei that has a bluish iris, these subspecies 
differ only in juvenile colour pattern, diag¬ 
nostic features lost during aging. Mauremys 
leprosa vanmeerhaghei (Oued Draa) is the 
only subspecies which should be easily rec¬ 
ognized as an adult. However, we encoun¬ 
tered that a pale blue iris may also occur in 
the Oued Noun (Fig. 4), and Keller & Busack 
(2001) highlighted additional inconsisten¬ 
cies regarding subspecific delineation. 

Our mtDNA sequences correspond to 
three of the newly described taxa (M. 1. ma- 
rokkensis, M. I. saharica, M. 1. vanmeer¬ 
haghei) and the nominotypical subspecies 
M. 1. leprosa (Schweigger, 1812) (type local¬ 
ity: South Spain). The subspecific allocation 
of our sample from Ceuta is unclear. Neither 


Schleich (1996) nor Bour & Maran (1999) 
attributed stripe-necked terrapins from that 
area to a specific subspecies. While it is 
beyond the scope of the present paper to 
discuss validity of any taxon, we wish to 
point out that we found the greatest sequence 
divergence in populations south of the Atlas 
Mts. This suggests some taxonomic differen¬ 
tiation. On the other hand, weak divergence 
of sequences from Ceuta, northern Morocco, 
and southern Spain questions the validity of 
subspecies described from northern Moroc¬ 
co. Likewise, we found south of the Atlas 
Mts. only weak differentiation between se¬ 
quences of turtles inhabiting the catchment 
areas of the Oued Draa (M. 1. vanmeerhaghei) 
and the Oued Noun (M. 1. saharica). In one 
sample from the Oued Dades, a tributary of 
the Oued Draa, we detected the same haplo- 
type as in our specimens from the Oued 
Noun, arguing for gene flow or incomplete 
lineage sorting and suggesting very recent or 
incomplete separation of both populations. 
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Taking the occurrence of these populations 
in isolated, in part temporary water bodies 
into account (Ducotterd & Bour 2002), the 
obvious differences in coloration and pattern 
of juvenile M. 1. saharica and M. 1. vanmeer- 
haghei (Fig. 5) could be caused by recent 
genetic bottlenecks. 

Our findings imply that while taxonomic 
differentiation within M. leprosa exists, it is 
currently overestimated with eight recog¬ 
nized subspecies, seven of which occur in 
Morocco. A complete phylogeography of M. 
leprosa will shed new light on its taxonomy 
and lead to a better understanding of the 
species’ biogeography. To decide whether 
vicariance, dispersal or both played a role in 
shaping current distribution patterns and in¬ 
fraspecific structure, we need additional sam¬ 
pling from North Africa and the Iberian pen¬ 
insula. 
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